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LONG-TERM  GOALS 

To  develop  formulas  to  evaluate  the  hydrodynamic  forces  acting  on  submerged  solid  bodies  using  only 
the  far  field  behavior  of  the  flow  past  the  body.  These  formulas  will  eliminate  the  need  to  resolve  the 
local  behavior  of  the  flow  in  the  vicinity  of  the  bodies.  Once  these  formulas  are  developed,  the 
interactions  between  submerged  solid  bodies  both  with  and  without  a  free  surface  may  be  analyzed 
without  the  need  to  resolve  the  local  fluid  velocity. 

OBJECTIVES 

Lighthill  [LI]  [L2]  used  a  decomposition  of  the  fluid  velocity  based  on  a  Newtonian  potential  to  relate 
the  Morison  equation  to  the  wave  induced  loads  on  submerged  bodies.  In  order  to  obtain  a  more  direct 
relation  between  the  hydrodynamic  forces  and  the  state  of  the  fluid,  different  representations  of  the  fluid 
velocity  are  required.  Accordingly,  our  objectives  are: 

•  To  improve  the  Lighthill  approach  by  considering  potentials  that  are  more  appropriate  for  studying 
fluid  flows,  such  as  the  Stokes  and  the  Oseen  fundamental  tensors. 

•  To  obtain  exact  or  asymptotic  formulas  that  relate  the  far  field  behavior  of  the  fluid  flow  to  the 
stress  vector. 

•  To  use  these  relations  for  the  evaluation  of  the  hydrodynamic  forces  on  submerged  solid  bodies. 
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APPROACH 


The  basic  principles  of  the  mechanics  of  fluids  and  solids  used  by  Lighthill,  [LI]  [L2]  were: 
a)  Decomposing  the  displacement  of  the  fluid  into  a  pure  straining  component  and  a  pure  rotational 
component;  and  b)  Obtaining  a  Helmholtz  decomposition  of  the  velocity  field  into  an  irrotational  and  a 
solenoidal  component.  The  Lighthill  decomposition  may  be  used  to  express  symbolically  the  loads  on 
submerged  bodies;  but  the  numerical  computation  of  these  loads  is  not  possible  due  to  absence  of 
boundary  conditions  for  this  decomposition.  One  of  the  limitations  is  the  Newtonian  potential.  This 
problem  may  be  avoided  in  the  linear  case  by  using  the  Stokes  or  Oseen  fundamental  solution  and  by 
focusing  instead  on  the  stress  vector. 

The  approach  will  develop  formulas  for  the  direct  evaluation  of  the  viscous  shear  stresses  for  steady 
nonlinear  flows.  A  solid  body  moving  with  a  constant  velocity  -U  in  a  fluid  at  rest  at  infinity  is 
analyzed.  The  Navier  Stokes  equations  for  the  steady  flow  of  an  incompressible  fluid  of  unit  density  are 
the  governing  equations  of  motion  for  the  fluid  occupying  the  region  exterior  to  the  solid  body  are 
given  by  [B]: 


[U  +  v)  •  Vv  =-pV  x  G)-  Vp  -gk  ;  V  •  v  =  0 


where  p  is  the  dynamic  fluid  viscosity,  «  is  the  vorticity,  and  p  is  the  dynamic  pressure.  On  the  surface 
of  the  solid  body  B,  the  no  slip  boundary  condition  is  assumed  to  hold  v(x)  =  -V  for  xc  B. 

The  approach  uses  potential  theory  identities  to  relate  the  far  field  behavior  of  the  dynamic  pressure  to 
the  total  forces  acting  on  the  submerged  solid.  This  is  a  departure  from  Lighthill’ s  approach  that 
considers  the  velocity  field. 

WORK  COMPLETED 

On  a  solid  body  with  angular  velocity  ,  the  stress  vector  may  be  written  as 

Tn  =  p(co  -  2  )xn  -  pn  (1) 


where  n  is  the  inward  directed  unit  vector  normal  to  the  surface  of  the  solid  body  B.  The  total  surface 
forces  acting  on  the  submerged  solid  body  are  given  by 
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For  the  Stokes  and  Oseen  problem,  Schuster  [S]  numerically  evaluated  the  stress  vector  on  the 
boundary  of  the  solid,  given  the  far  field  velocity  by  solving  a  Fredholm  equation  of  the  first  kind.  The 
kernel  of  the  integral  equation  is  obtained  from  the  fundamental  solution  tensor  of  the  corresponding 
equation.  Relations  between  the  velocity  at  infinity  and  the  stress  vector  appear  to  be  limited  to  only  the 
linear  problems.  Consequently,  the  dynamic  pressure  will  be  considered. 


For  the  nonlinear  steady  flow  problem,  the  following  new  formula  for  dynamic  pressure  was  derived: 


4  p(x)  =  Jr(y)n(y)-V^-L_rfS(y)  -  J(v - V) v - V^— L_rfy 
-  r(U'V)vn(y)  — - — dS(y)  -  fg-^-dS( y) 

{  x~y  {  x~y 


(2) 


The  relation  between  the  far  field  behavior  of  the  dynamic  pressure  and  the  stress  vector  may  be 
obtained  as  follows.  Multiplying  (2)  for  p  by  x/R  and  integrating  over  the  sphere  of  radius  R,  the  last 
two  integrals  vanish.  In  addition 


/ 


x  =R 


— ®V- 
R 


dS(x)  =  -4  1  +  0 
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where  I  is  the  3x3  identity  matrix.  The  volume  integral  in  (2)  for  the  pressure  formula  decays  according 
to  the  results  given  by  Finn  [F].  The  hydrodynamic  forces  on  the  submerged  solid  body  may  then  be 
computed  from 
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The  details  of  this  derivation  are  given  in  [GHT]. 

RESULTS 

The  numerical  calculations  by  Schuster  [S]  for  the  linear  case  show  that  it  is  possible  to  relate  the  solid 
body  velocity  to  the  stress  vector.  The  relations  that  one  obtains  require  the  fundamental  solution  of  the 
Stokes  or  Oseen  problem,  and  may  not  be  accomplished  using  a  Newtonian  potential  as  suggested  by 
Lighthill  [LI]  [L2].  The  numerical  solution  of  the  integral  equations  that  arise  requires  the 
characterization  of  the  null  space  that  is  given  by  Fisher,  Hsiao  and  Wendland  [F],  [HW],[FHW]. 

In  the  nonlinear  regime,  it  does  not  appear  possible  to  relate  the  solid  body  velocity  to  the  stress  vector. 
In  part,  this  is  due  to  the  coupling  of  the  kinematic  boundary  conditions  for  the  scalar  and  vector 
potentials  used  in  the  Helmholtz  decomposition  of  velocity  vector  fields.  Consequently,  relations 
between  the  stress  vector  and  thermodynamic  variables  ( e.g .,  dynamic  pressure)  are  required.  An 
example  of  this  relation  for  steady  flow  is  given  by  (3). 

IMPACT/APPLICATIONS 

The  new  formulas  (2)  and  (3)  appear  to  provide  a  method  for  computing  fluid  drag  forces  that  is  not 
possible  using  the  Lighthill  velocity  decomposition. 


The  significance  of  (1)  is  that  the  stress  vector  may  be  related  pointwise  to  the  vorticity  of  the  fluid. 
Consequently,  (1)  may  be  used  to  derive  rigorous  or  semi-rigorous  models  for  the  study  of  sediment 
transport. 

TRANSITIONS 

These  results  lead  to  the  development  of  efficient  numerical  methods  for  solving  Fredholm  integral 
equations  of  the  first  kind.  Equation  (3)  offers  an  appealing  alternative  to  far  field  velocity  [LI1]  or 
vorticity  [R]  methods  for  computing  the  hydrodynamic  force  on  submerged  solid  bodies  that  could  have 
substantial  impact  on  the  analysis  of  loading  on  structures. 

RELATED  PROJECTS 

The  results  of  this  research  could  be  adapted  and  reformulated  to  analyze  sediment  transport.  These 
results  could  also  be  used  to  determine  the  boundary  data  needed  in  numerical  algorithms  based  on 
vortex  elements  for  computational  fluid  mechanics. 

REFERENCES 

[B]  G.  Batchelor,  1973:  An  Introduction  to  Fluid  Dynamics ,  Cambridge  University  Press. 

[F]  T.  Fischer,  1987:  “A  boundary  integral  method  for  the  numerical  computation  of  the  forces  exerted 
on  a  sphere  in  viscous  incompressible  flows  near  a  plane  wall,”  J.  of  Applied  Mathematics  and  Physics 
(ZAMP),  Vol  38.  pp  339-365. 

[GHT]  R.  Guenther,  R.  Hudspeth,  E.  Thomann,  “The  evaluation  of  the  hydrodynamic  forces  on  a  solid 
body  from  the  far  field  pressure  field,”  In  preparation. 

[HW]  G.  Hsiao  and  W.  Wendland,  1977:  “A  Finite  element  method  of  some  integral  equations  of  the 
first  kind,”  J  of  Mathematical  Analysis  and  Applications ,  Vol.  58.  pp  449-481. 

[FHW]  T.  Fischer,  G.  Hsiao,  and  W.  Wendland,  1985:  “Singular  perturbations  for  the  exterior  three 
dimensional  slow  viscous  flow  problem,”  J  of  Mathematical  Analysis  and  Applications,  Vol.  110.  pp 
583-603. 

[FH]  M.  Fonguet  Higgins,  1992:  “Capillary  rollers  and  bores,”  J.  of  Fluid  Mechanics,  Vol.  240  pp 
659-679. 

[FI]  J.  Fighthill,  1986:  “Fundamentals  Concerning  Wave  Foading  on  Offshore  Structures,”  J.  Fluid 
Mechanics,  Vol.  173.  pp  667-681. 

[F2]  J.  Fighthill,  1986:  An  Informal  Introduction  to  Theoretical  Fluid  Dynamics,  Oxford  University 
Press. 

[R]  D.  Rockwell,  Ocean  Engineering  and  Marine  Systems  1997  Program,  ONR,  pp.  54-63. 

[S]  M.  Schuster,  1988:  “Computation  of  the  Stresses  on  a  Rigid  Body  in  Exterior  Stokes  and  Oseen 
Flows,”  MS  Thesis,  Dept.  Mathematics,  Oregon  State  University,  June. 

PUBLICATIONS 

Wang,  H.  and  Guenther,  R.B.  1997:  “Calculation  of  the  Far  Field  Finite  Depth  Green's  Function,”  J. 
Ocean  Engineering,  24  (1),  83-94. 

Foster,  D.F.,  Guenther  R.B.,  Holman,  R.A.  1998:  “An  Analytic  Solution  for  the  Wave  Bottom 


Boundary  Layer  Governing  Equation  Under  Arbitrary  Wave  Forcing,”  to  appear  in  J. Ocean  Engr. 


Guenther,  R.B.  and  Lee  J.W.  1998:  “Heat  Conduction  with  Radiating  Boundary  Conditions,”  Journal 
of  Computational  and  Applied  Mathematics,  88,  119  -  124. 

Guenther,  R.B.  and  Lee,  J.W.  “Boundary  Value  Problems  for  a  Class  of  Integro-differential  Equations 
and  Applications,”  to  appear  in  Journal  of  Computational  and  Applied  Mathematics. 

Bowline,  C.,  Hudspeth,  R.T.,  Guenther,  R.B.,  “Are  Cross  Waves  Chaotic?”  to  appear  in  J.  Applicable 
Analysis. 

Hudspeth,  R.T.,  Borgman,  L.E.,  and  Samorsom,  B.,  '’Conditional  Simulation  of  Laboratory  Waves,”  to 
appear  in  Journal  of  Waterways,  Port,  Coasted  and  Ocean  Engineering. 

Tuah,  Hang,  Chen,  M-C,  and  Hudspeth,  R.T.  1998:  “Discussion  of  Numerically  Simulating  Non- 
Gaussian  Sea  Surfaces”  by  B.  Vanhoff,  S.  Elgar  and  R.T.  Guza,  Journal  of  Waterways,  Port,  Coastal 
and  Ocean  Engineering,  Vol.  124  (6),  Nov/Dec,  pp.  335-336. 

Hudspeth,  R.T.,  “Discussion  of  an  Improvement  to  Stokes  Nonlinear  Theory  for  Steady  Waves” 
by  V.  Karambas  and  C.G.  Routitas  to  appear  in  Journal  ofWatem’ays,  Port,  Coasted  einel  Oceein 
Engineering. 


